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Discrete optical solitons were first suggested theoretically as stationary nonlinear 
localized modes of a periodic array of weakly coupled optical waveguides. Because 
the use of discrete solitons promises an efficient way to control multi-port nonlinear 
switching in a system of many coupled waveguides, this field has been extensively 
explored theoreticallypEEEl More importantly, the discrete solitons have also been 
generated experimentally in fabricated periodic waveguide structures)^ 

The most common theoretical approach to study discrete optical solitons in 
waveguide arrays is based on the decomposition of the electric field in the periodic 
structure into a sum of weakly coupled fundamental modes excited in each waveguide 
of the array. According to this approach, the wave dynamics can be described by 
an effective discrete nonlinear Schrodinger (DNLS) equation, that possesses spatially 
localized stationary solutions in the form of localized modes of a lattice model. Many 
properties of the discrete optical solitons can be analyzed in the framework of the 
DNLS equation.™ll 

One of the major problems for achieving the controllable multi-port steering of 
discrete optical solitons in waveguide arrays is the existence of an effective periodic 
potential which appears due to the lattice discreteness, known as the Peierls-Nabarro 
(PN) potential. It represents the energy cost associated with a shift of a nonlinear 
localized mode by a half of the waveguide spacing.- Its magnitude can be roughly 
estimated as ~ |^4| 4 , where A is the soliton amplitude. As a consequence of this po- 
tential, a narrow large-amplitude discrete soliton does not propagate in the lattice and 
it becomes trapped by the array. Several ideas to exploit the discreteness properties of 
the array for the optical switching were suggested,'*^' including the demonstration of 

2 



the output channel selection for the multi-port devices. However, the soliton steering 
and switching is well controlled only in the limit of broad beams whereas the soli- 
ton dynamics in highly discrete arrays has been shown to be more complicated.^ In 
this Letter, we suggest a "discreteness engineering" approach and demonstrate how 
to achieve highly controllable multi-port soliton switching in the arrays by a desired 
integer number of waveguides, the so-called "digital soliton switching" . 

We consider a standard model of the waveguide arrays with a modulated coupling 
described by the normalized DNLS equation of the form, 

du 

i ~T L + V n+lUn+l + + 7l«n| Vi = 0, (1) 

dz 

where u n is the effective envelope of the electric field in the n-th waveguide and z is 
the propagation distance. Unlike the standard modelsf^EEl the coupling V n between 
two neighboring guides is assumed to vary, either through the effective propagation 
constant or by changing the spacing between neighboring guides. The parameter 
7 = uj n 2 / (cA ef [) is the effective waveguide nonlinearity associated with the Kerr 
nonlinearity of the core material. 

The steering and trapping of discrete solitons in the framework of the model 
have been analyzed in many studies. Being kicked by an external force, the dis- 
crete soliton can propagate through the lattice for some distance, but then it gets 
trapped due to the effect of discreteness. For a larger force, the output soliton posi- 
tion fluctuates between two (or more) neighboring waveguides making the switching 
uncontrollable.^ Here, we suggest to modulate the waveguide coupling in order to 
achieve a controllable output and to engineer the switching results. The key idea is 
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to break a symmetry between the beam motion to the right and left at the moment 
of trapping; this allows the elimination of chaotic trapping observed in homogeneous 
arrays.^ In this way, we achieve a controllable digitized switching where the contin- 
uous change of the input beam amplitude results in a quantized displacement of the 
output beam by an integer number of waveguides. 

We have tested different types of modulation in the array parameters and the 
corresponding structures of array super-lattices. An example of one of the optimized 
structures is shown in Fig. 1, where we modulate the coupling parameter V n in a step- 
like manner. We also notice that the use of a linear ramp potential (e.g. of the form 
V n = an) for this purpose does not lead to an effective switching but, instead, makes 
the soliton switching even more chaotic due to the phenomenon of Bloch oscillations 
which become randomized in the nonlinear regime.^ 

We select the input profile in the form of a narrow sech-like beam localized only 
on a few waveguides, 

u n (0) = A sech[A{n - n c )/V2] e ~ lk{n -^\ (2) 

for n — n c = 0, ±1, and u n = 0, otherwise. For the particular results presented below, 
we select the array of 101 waveguides and place the beam at the middle position, 
n c = 50. The maximum normalized propagation distance used in our simulations is 
z max = 45 (in units of the coupling length). 

Parameter k in the ansatz (j2J) has the meaning of the transverse steering velocity 
of the beam, in analogy with the continuous approximation. It describes the value 
of an effective kick of the beam in the transverse direction at the input, in order 
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to achieve the beam motion and shift into one of the neighboring (or other desired) 
waveguide outputs. 

In our simulations, we control the numerical accuracy by monitoring two con- 
served quantities of the model (J2J), the soliton power P = J2 n \ u n(z)\ 2 , and the system 



The input condition (J2J) is not an exact stationary solution of the discrete equation 
(JIJ even for k = 0, and as the input kick (k ^ 0) forces the soliton to move to the right 
(k < 0) or left (k > 0), the motion is accompanied by some radiation. The effective 
lattice discreteness can be attributed to an effective periodic PN potential. Due to 
both the strong radiation emission and the PN barrier which should be overtaken in 
order to move the beam, the discrete soliton gets trapped at one of the waveguides 
in the array, as shown in Fig. 2. In most cases, the shift of the beam position to the 
neighboring waveguide is easy to achieve, as shown in many studies.^ However, the 
soliton switching becomes rather complicated and even chaotic when the kicking force 
becomes stronger. 

We have studied many different regimes of the soliton multi-port switching in the 
array and revealed that the most effective switching in a desired waveguide position 
(i.e. desired output) can be achieved by varying the coupling between waveguides, 
either through the effective propagation constant or by changing the spacing between 
neighboring guides, as shown in Fig. 1. This V n profile was obtained after performing 
a numerical sweep in V n and A for fixed momentum k. In this case, the selection 
of a finite value of the steering parameter k allows to switch the whole beam into a 
neighboring waveguide, as shown in Fig. 2, with only a small amount of radiation. By 
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decreasing the amplitude of the input pulse at a fixed value of the steering parameter, 
fixed to be say k = ±0.9, it is possible to achieve self-trapping of the soliton beam 
at some (short) distance from the initial center at different waveguide position. Due 
to the step-like modulated coupling, we create a selection between the beam motion 
to the right and left at the moment of trapping thus suppressing or eliminating the 
chaotic trapping observed in homogeneous waveguide arrays. In this way, we achieve 
a controllable digitized nonlinear switching where the continuous change of the am- 
plitude of the input beam results in a quantized displacement of the output beam by 
an integer number of waveguides. Consequently, for the parameters discussed above 
we observe almost undistorted switching up to eleven waveguides, and Fig. 3 shows 
an example of the digital soliton switching to the eleventh waveguide.^ 

Figure 4 gives a summary of the results for the parameters discussed above; it 
shows the discrete position of the soliton at the output as a function of the input beam 
amplitude, for two fixed values of the steering parameter k = ±0.9. In a remarkable 
contrast with other studies, the coupling modulation allows a controllable digitized 
switching in the array with very little or no distortion. The figure also shows a slight 
asymmetry in the final displacement, depending on whether the beam is kicked uphill 
or downhill. 

If we were to use five guides instead of three for the input beam, one could 
expect a smaller amount of radiation emitted. However, this would imply a longer 
distance before the beam gets trapped by one of the waveguides in the array due to 
the effective Peierls-Nabarro potential. Also, this means that one could in principle 
switch the soliton beam to any desired waveguide in the waveguide array, no matter 
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how far; it would be just a matter of choosing an initial beam wide enough, i.e., closer 
to the continuum limit (in addition to the optimization of the coupling in a step-wise 
manner) by removing the random selection between the directions and suppressing 
the beam random switching. 

Another observation is that the sech-like initial profile is not really fundamen- 
tal. A (kicked) nonlinear impurity-like profile of the form u n (0) = A[(l — A 2 )/(l + 
^2^j|n-n c |/2 e xp[— ik(n — n c )] will also show similar behavior, as our additional compu- 
tations show. The reason for this behavior seems to rest on the observation that, for 
any system with local nonlinearity, a narrow initial profile will effectively render the 
system into a linear one containing a small nonlinear cluster (or even a single site); 
the bound state will therefore strongly resemble the one corresponding to a nonlinear 
impurity.^ 

In conclusion, we have suggested a novel approach to achieve a digitized switching 
in waveguide arrays by using the concept of discrete optical solitons. Our approach 
involves a weak step-like modulation of the coupling strength (or, equivalently, dis- 
tance between the waveguides) in the arrays with the period larger than the waveguide 
spacing. Such a super-lattice structure allows the modification of trapping properties 
of the waveguide array due to discreteness which in turn permits engineering of the 
strength of the effective trapping potential. We have demonstrated numerically the 
controlled switching up to eleven waveguides in the arrays by using very narrow input 
beams localized on three waveguides only. 
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Figure Captions 



Fig. 1. Example of the optimized modulation of the propagation constant V n in 
the waveguide array. 

Fig. 2. One-site switching of a discrete soliton in the waveguide array with the 
modulated coupling shown in Fig. 1. 

Fig. 3. Discrete switching by eleven sites in the waveguide array modulated ac- 
cording to Fig. 1. 

Fig. 4. Soliton switching in a waveguide array with an optimized coupling. Shown 
is the soliton output displacement as a function of the input beam amplitude (A). A 
step size of 0.0005 separates consecutive points. 
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Fig. 1. Example of the optimized modulation of the propagation constant V n 
in the waveguide array. 
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Fig. 2. One-site switching of a discrete soliton in the waveguide array with the 
modulated coupling shown in Fig. 1. 
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Fig. 3. Discrete switching by eleven sites in the waveguide array modulated 
according to Fig. 1. 
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Fig. 4. Soliton switching in a waveguide array with an optimized coupling. 
Shown is the soliton output displacement as a function of the input beam 
amplitude A. 
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